Schmidt-like coherent mode decomposition and spatial intensity correlations of 

thermal light 
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We experimentally study the properties of coherent mode decomposition for intensity correla- 
tion function of quasi-thermal light. We use the technique of spatial mode selection developed for 
studying transverse entanglement of photon pairs, and show that it can be extended to characterize 
classical spatial correlations. Our results demonstrate the existence of a unique for a given thermal 
source basis of coherent modes, correlated in a way much resembling the Schmidt modes of spatially 
entangled photons. 



I. INTRODUCTION 

Entanglement in spatial degrees of freedom of quan- 
tum light, such as pairs of photons generated in sponta- 
neous parametric down-conversion (SPDC), is currently 
an object of active research. These studies are motivated 
by applications in quantum information science, where 
infinite-dimensional Hilbert space of spatial states of pho- 
tons offers attractive capabilities for high-dimensional 
quantum state engineering. The effective dimension- 
ality of accessible state space is closely related to de- 
gree of spatial entanglement for bipartite states, mak- 
ing its experimental quantification an important task. 
One of the most successful approaches makes use of a 
Schmidt decomposition. As stated by Law and Eberly in 
their seminal work a unique basis of coherent modes 
{\ u k) , \vk)} exist, such that a spatial state of biphoton 
pairs 1^12) = / dxidx2^(xi, a; 2 )a^ (xi)aJ(£ 2 ) |vac) can 
be decomposed as |\E r ) = J2k yf^k\ u k) \ v k), with these, 
so-called, Schmidt modes being eigenvectors of reduced 
single-particle density matrices and - corresponding 
eigenvalues. A remarkable feature of this decomposition 
is its single-sum form, implying perfect one-to-one corre- 
lations between Schmidt modes. These correlations were 
studied in several recent works , and number of sig- 
nificant terms in the decomposition is routinely used as 
an entanglement quantifier 0, [EMT^ . 

A natural question is whether similar approach may 
be used to study spatial correlations of other origin, for 
example, correlations in spatially multi-mode classical 
light? The form of Schmidt decomposition, and the way 
explicit expressions for Schmidt modes are derived in the 
case of SPDC biphotons [J [l3j resembles the well known 
notion of a coherent mode decomposition for spatial cor- 
relation functions in classical statistical optics p^4l6| . 
Here we show, that this similarity does not only reside 
in the form of mathematical equations, but in underly- 
ing physics as well. Namely, for a given source of light 
with quasi-thermal statistics there exist a unique phys- 
ically distinguished set of coherent modes, correlated in 
much the same way as Schmidt modes of SPDC. Similar 



experimental settings would provide similar correlations, 
with only a quantitative difference in visibility. The situ- 
ation here resembles that with ghost imaging, which was 
for some time believed to be possible only due to en- 
tanglement [i~7j and, thus only with quantum light. It 
is however well understood now, that spatial correlations 
required to obtain ghost images should not necessarily be 
of quantum origin, and a purely classical thermal light 
may be used as well [HI, Il9j . In the present work we 
try to push this analogy further and show, that a frame- 
work of Schmidt decomposition may also be generalized 
to classical states of light. 

The paper is organized as follows: in section II we 
review some properties of coherent mode decomposition 
for a Gaussian Shell-model of quasi-thermal light and 
derive the expressions for measured correlation functions, 
experimental realization is described in section III, and 
we end with discussion of results and conclusion in section 
IV. 



II. COHERENT MODE DECOMPOSITION FOR 
QUASI-THERMAL LIGHT 

Spatial coherence properties of classical light are 
described by coherence function G^^ri,^) = 
(E*(ri)E(r2)), where we neglect the polarization 
and use scalar complex amplitudes E(r) to describe 
electromagnetic field. One of the results of classical 
coherence theory states that it can be decomposed in 
series of coherent mode functions as 



G( 1 )(r 1 ,r 2 ) = ^A„C(ri)^(r 2 ), 



(1) 
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where <fi n (r) are eigenfunctions of an integral operator 
with a kernel G*- 1 ) (ri , r 2 ) and A„ are corresponding eigen- 
values [l6j . 

Higher order correlation functions, such as the inten- 
sity correlation function G^ 2 '(ri, r 2 ) = (I(ri)/(r 2 )), may 
be accessed via Hanbury Brown- Twiss (HBT) interfer- 
ometry, which is exactly the experimental setting com- 
monly used for studying spatial entanglement of pho- 
ton pairs. In the biphoton case the fourth-order cor- 
relation function G^ 2 - I (ri,r 2 ) is measured by counting 
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coincidences between photocounts of two single-photon 
detectors positioned at points and r 2 , respectively: 
R c oc (at(n)at(r 2 )a(ri)a(r 2 )) = |*(n,r 2 )| 2 , with a(r) 
being the photon annihilation operator and 4'(ri,r 2 ) - 
the biphoton amplitude. A direct analogue of coherent 
mode decomposition (p} for the biphoton amplitude is 
exactly the Schmidt decomposition. To see whether this 
analogy goes beyond simply similar mathematical expres- 
sions, let us consider the behavior of classical light field 
in the HBT scheme with spatial mode filters in the arms 
of the interferometer. Such filters are usually composed 
of phase holograms followed by single mode fibers and 
bucket detectors as shown in Fig. [T] The main feature 
of Schmidt decomposition - its single sum form, mani- 
fests itself in absence of photocounts in such setup when 
appropriate mode filters are used to select orthogonal 
Schmidt modes in different arms of the setup 0. We 
will show below, that similar effects may be observed 
with classical quasi-thermal light. 

For quasi-thermal light the intensity correlation func- 
tion G< 2 )(ri , r 2 ) can be expressed in terms of the second- 
order correlation function G^^(ri,r 2 ) using Siegert rela- 
tion: 



G( 2 >(n,r 2 ) = (/(n))(/(r 2 )) + |G (1) (r!,r 2 )p 



(2) 



If the spatial mode filters in the two arms of the HBT 
scheme are described by impulse response functions 
hi(ri, r'j) and /i 2 (r 2 , r 2 ), where coordinates r corresponds 
to the source plane and r' - to the detection plane, the 
transformed intensity correlation function takes the form 

MM- 

G( 2 )(r' 1 ,r 2 ) = (7(r' 1 ))(7(r 2 )) 
+ |/dr 1 /dr 2 ^(r 1 ,r' 1 )/ l2 (r 2 ,r 2 ) G«(n,r 2 )| 2 , (3) 

where I(r' 12 ) = \ J dr^ 2 /ii i2 (r^ 2 , r[ 2 )i?(ri, 2 )| 2 is the in- 
tensity distribution in the detection plane. One can al- 
ways choose the mode filters to satisfy the relation: 

(r x ) oc e w t 6 



dr 1 /i 2 fe) (r 2 ,r 2 )^„(r 2 ) 



"f Sr, 



(4) 



where Wf is the waist of the fundamental Gaussian mode 
of the fiber. This choice of propagators corresponds to 
"projection" on modes <fi m (ri) and </>fc(r 2 ). Substituting 
(H|) and ([3]) into © and taking into account that bucket 
detectors placed after the single mode fiber integrate the 
signal over the whole detection plane, we obtain the fol- 
lowing expression for coincidence to single counts ratio: 



R 



p (m) p (fe) 



(2) 
9 (jn,k) 



OC 1 + <L 



(5) 



where 



(2) 
9(m,k) 



fdri(/(rO> fdr 2 (/(r 2 )) 



is the normalized second order correlation function. It is 
clear from (|5|), that spatial correlations between appro- 
priately chosen coherent modes of thermal light show the 
same correlation features as Schmidt modes of entangled 
photons - namely, they are only pairwise correlated. The 
difference turns out to be rather quantitative, than qual- 
itative, similar to the situation with thermal and quan- 
tum ghost imaging - the achievable visibility defined as 



V = (g 



(2) 
(m,m) 



S(m,n))/(ff 



(2) 

(m,m) 



(2) 

(m,n) 



), Vm ^ n can 



not exceed 1/3 for the thermal light, while it can reach 
values arbitrarily close to unity for biphotons. 

Let us now derive the explicit expressions for coher- 
ent modes constituting the decomposition (TTJ) for a sim- 
ple model of partially spatially coherent thermal light. 
Following [l4[ we describe light from the quasi-thermal 
source with a Gaussian Shell-model for the second order 
coherence function: 



G (1) (xi, x 2 ) = y / I{xi,ui)I(x2,uj)n(xi - x 2 ,uj), (6) 

where spectral intensity distribution and degree of spatial 
coherence arc taken in the form 



I{x,w) = A(w)exp[-x 2 /2cr|(w)], 



(7) 



fi(xi - X2,w) = cxp[-(xi - x 2 ) 2 /2cr 2 (w)], (8) 

with X\, x 2 being the transverse coordinates, A(uj) a pos- 
itive normalization constant, cri(uS) - a beam waste and 
cr^ui) - its coherence radius. For the sake of simplicity 
we will first consider a one-dimensional source, gener- 
alization to an experimentally relevant two-dimensional 
case is straightforward. We will also assume the light to 
be monochromatic and omit the frequency dependence. 

In the case of partially spatial coherent light described 
by a Gaussian Shell-model, one can find the eigenfunc- 
tions in a closed analytical form [l4|, [H| : 



1/4 



\/2 n nl 



i?„(xV2c)exp(-cx 2 ), (9) 



A„ = A 



a + b + c 



1/2 



(10) 



where 



r, b 



i 7 , c= y/a 2 + 2ab, (11) 



4a 2 ' 4a 2 

and H n (x) are Hermite polynomials. Eigenvalues in this 
coherent mode decomposition are determined by the ra- 
tio of coherence radius to beam waist B = — : 



1 



~ ~~ V/? 2 /2 + l + /?[(/?/2) 2 + l]V; 



(12) 



In this simple, but experimentally relevant model the 
cigenmodes turn out to be familiar Hermite- Gaussian 
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(HG) modes, which were also found to be Schmidt modes 
of SPDC biphotons under a double-Gaussian approxima- 
tion for the biphoton amplitude [HEl, EH- We note, that 
similar decomposition may be carried out in polar coordi- 
nates giving rise to coherent modes with orbital angular 
momentum (OAM) plM^. Spectrum of OAM eigenval- 
ues in this decomposition was measured experimentally 
in (2~H with interferometric technique. 

III. EXPERIMENTAL REALIZATION 

The experimental setup is shown in FigUJ We used 
a well-known technique to create a pseudothermal light 
source with controlled spatial coherence (2f|. In this 
scheme a beam of a He-Nc laser with wavelength A = 
632.8 nm is scattered on a slowly rotating groundcd-glass 
disc (GG). It is important to make the characteristic time 
of variation of the scattering pattern larger than the de- 
tection time Tdet ~ 2ns, defined by the coincidence cir- 
cuit time window. The disc was located at the focus 
of an effective lens consisting of two microscope objec- 
tives 01 and 02, forming a collimated beam. Spatial 
coherence properties of the beam are controlled by iris 
apertures S and S2, determining the coherence radius <7 M 
and the beam waist 07, respectively. The detection part 
of the setup is an HBT scheme, consisting of a 50/50 
non-polarizing beam-splitter; spatial mode filters - a re- 
flective phase-only LCoS SLM in the transmitted channel 
and a glass step-mask in the reflected one, followed by 



single-mode fibers SMF\_2 placed in the focal planes of 
8x microscope objectives Ci^; electrical signals from two 
single-photon detectors Di^ (Perkin-Elmer) were fed to 
a coincidence curcuit. The model of SLM used was lim- 
ited to 0.87T phase shifts only, forcing us to use double- 
reflection as shown in Fig. Q] Polarizers P2,3 were used 
to specify the correct polarization required for operation 
of polarization sensitive SLM matrix and to remove the 
undiffractcd components, while polarizer P\ was used to 
control the overall intensity. 

He-Ne laser 




FIG. 1. Experimental setup. (See text for detailed de- 
scription) 

The propagator for each arm of the HBT scheme de- 
scribed above is 



h(r, r') oc exp 



(r') 2 



"'7 



dr" cxp(— i - ) x exp iarg ^H m (v / 2ca;)H n (v / 2c2/)J 



x exp 



(r" - r f f 



(13) 



where k = 2tt/X is the wave number, / - the focal length 
of the microscope objectives Ci,2, Wf - the waist of the 
fundamental gaussian mode of the fiber, - the fiber 
tip displacement in the focal plane. This propagator de- 
fines the detection modes, which for our purposes should 
be matched with eigenmodes given by Q. The mode 
matching is accomplished by choosing the focal length 

f = \ -p- 1 in this case the detection mode waist coin- 

J V 4c ' 

cides with the waists of the eigenmodes: 



k 2 w 2 j 
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The mode filters were tested with a coherent single 
mode laser beam. For this purpose the grounded-glass 
disc was replaced by a mirror and radiation of a 650 
nm diode laser was coupled into SMF\ fiber in "back- 
ward" direction (see FigfT]). The after passing through 
both arms of the setup the beam was detected in the 
reflected arm. One may consider this procedure to be 



the direct implementation of Klyshko's " advanced- wave 
picture" (26|. The intensity at the output of the fiber 
is proportional to the convolution of HG m „ and HGrjo 
modes 



HG^OHGoofo-Ode 



(14) 



where Xf was the horizontal coordinate of the tip of 
SMF2. Convolutions for the first two modes HGrjo and 
HG10 are shown in Fig. [3] Data shown in the Figure 
[2ja) were used to determine the width of the fiber mode 
(wf = 3.69 ± 0.04 /Ltm). Phase masks used for mode 
conversion had step-like phase profiles shown in Fig. [3] 
Sizes and positions of the phase masks relative to the 
beam were determined experimentally following the op- 
erational criterion of minimizing the counting rate in the 
central position of the SMF2 fiber. 

k 2 w 2 ' 

Satisfying the condition -%jt- = c is essential to ob- 
tain the desired matching of the detected modes and 
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FIG. 2. Hermite-Gaussian mode filtering. Convolution 
of HGoo (a) and HGio (b) modes with gaussian fiber 
mode. Solid red lines are: Gaussian fit used to estimate 
detection mode waist (a), and HGio function with the 
same waist (b). 
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m=0, n=2 


H m=2 f n=2 











FIG. 3. Phase masks corresponding to Hermite-Gaussian 
modes of different orders. Masks parameters are deter- 
mined experimentally (see text for details). 



the eigenmodes of coherence function. In the actual ex- 
periment the detection mode waist was fixed, while the 
mode-matching was done by changing the parameters of 
the input beam. Width of the apertures S and S2 de- 
termine the coherence radius and the beam waist 07, 
correspondingly, and thus define the eigenmodes width 
c through (fTTj) . One can show by direct calculation 

that the value of normalized intensity coherence function 
(2) 

9(m o)( r / = 0) = 1 for even m only if the mode- matching 
condition is satisfied, otherwise it is always larger than 
1, as shown in Fig. @] This fact was used as an opera- 
tional criterion of good mode-matching. The correspond- 
ing beam waist and coherence radius were measured to 
be 01 = (2.3 ± 0.1) mm and 07, = (0.57 ± 0.02) mm, 
corresponding to the value of (3 = 0.24 ± 0.02. 

With the mode-matching conditions satisfied, the spa- 
tial mode filters in both arms of the setup act as projec- 
tors on eigenmodes of coherent mode decomposition (jTJ) . 
In this case the single counts rate of a detector in one 
of the arms is proportional to partial intensity of a con- 
stituent mode selected by the filter: 



(ro) 



oc 



rir' 



drh^(r,r',)E(v) 



oc X Tl 



(15) 



It provides a way to measure the eigenvalues in (TTJ exper- 
imentally. The distribution of eigenvalues, normalized to 
Aq is shown in Fig. [5] We find excellent agreement with 
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FIG. 4. Convolutions of the gaussian detection mode 
with HG20 modes of different waists, (a) The waist c is 
"incorrect" , the mode-matching condition is not satisfied, 
normalized fl^o)^/ = 0) > 1. (b) The waist is chosen 
appropriately, normalized intensity correlation function 
reaches the minimal value of unity. 



theoretical predictions obtained from (TT2"j) using the in- 
dependently measured value j3 = 0.24 ± 0.02. 



(b) 




o OA 

^ n 



mode index, m 

FIG. 5. Experimentally measured eigenvalues of coher- 
ent mode decomposition, (a) Two-dimensional distribu- 
tion of normalized eigenvalues A mn /An (indexes m and 
n correspond to HG m o and HGon modes, respectively), 
(b) One-dimensional section (n = 0) of the distribution. 
Blue bars are experimental data, black circles are values 
calculated from (|12|) using independently measured value 
f3 = 0.24 ± 0.02. Error bars correspond to error in mea- 
sured value of /3, error bars for experimental points are 
smaller than point markers. 

The spatial shape of modes © selected by correspond- 
ing mask in the transmitted channel is revealed in the 
dependence of intensity correlation function on fiber dis- 
placement in the reflected channel. Indeed, if a phase 
mask corresponding to m-th eigenmode is installed in the 
transmitted channel with no mask in the reflected, using 
([3]) one can find the following expression for a normalized 
intensity correlation function: 



|A m a m (r/)|' 



fc 



(16) 



where 17 is the position of the single mode fiber in the 
reflected channel, and 



«fc( r /) = / drr> fc (r)r> (r - 17). 



(17) 



Experimental data for HG m rj modes are shown in Fig. [B] 
It is clear, that the selected modes arc indeed orthogo- 
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(2) 

nal in the sense, that sLjq)^) = 1, Vm ^ 0, consis- 
tently with ([5]). Dependence on the fiber displacement 
is also in qualitative agreement with (|16[) . However for 
large displacements we observe anomalously high values 

( 2) 

of g^ j(r/) for large m. This is most probably an arte- 
fact of our method for mode selection based on using 
phase-only step masks, inevitably causing higher-order 
mode contributions and compromising mode purity for 
large displacements. 
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FIG. 6. Spatial structure of eigenmodes: dependence of 

(2) 

normalized gY^ j function on fiber displacement in the 
reflected channel. Unity values for zero displacement cor- 
respond to absence of correlations for modes of different 
order, predicted by theory. 

For zero displacement and low-order modes we find the 
detected modes to be close to Hcrmite- Gaussian eigen- 
modes, and thus they are pairwise correlated, as illus- 
trated by Fig. [7] We expect g^ n \{pf = 0) to behave 
in accordance with expression (0, i.e. to be essentially 
unity for modes of different order, and larger than unity 
only for modes of the same order. Experimental results 
of Fig. [7] confirm this statement, although some amount 
of anomalous correlations is observed for modes of high 
order. These correlations are more pronounced on Fig. [7] 
(b), which may be partially explained by much poorer 
quality of the glass phase-masks used in the reflected 
channel. Nevertheless, one can clearly observe the good 
agreement with expected correlation properties of quasi- 
Schmidt coherent modes. 



IV. DISCUSSION AND CONCLUSION 

The main source of discrepancies between our experi- 
mental results and theoretical expectations for ideal pro- 
jective measurements in Schmidt-like basis is the usage 
of non-ideal mode filters. Technical limitations of the 
SLM used do not allow to use it for amplitude mod- 
ulation, which is necessary for perfect transformations 
of Hcrmitc-Gaussian modes. Nevertheless, our results 





FIG. 7. Mode cross-correlations: 



1 values mea- 



( a ) 9(m,0) 

sured with HG m o modes selected in the transmitted arm and 
HGoo mode in the reflected one; (b) ^ — 1 values, cor- 
responding to HG m o modes filtered in the transmitted arm, 
and HGoi mode selected in the reflected arm. 



are very similar to what is expected for exact Hcrmitc- 
Gaussian filters. Let us address this issue quantitatively. 
One can define fidelity 

N,M 



12 V ^mn^ 



exp 
mn 



F = 



i 



M,JV 



(M, N N 
\ ex P 
/ i ^mn 



(18) 



as a measure of deviations of (normalized) distribution of 
experimentally obtained eigenvalues from the prod- 
uct of two (normalized) theoretical distributions (fT2|) for 
a two-dimensional Gaussian Shell-model of incoherent 
light. The value of this quantity for our data is depen- 
dent on the upper limits {M, N}, i.e. on the maximal 
order of modes, that are taken into account. This de- 
pendence is shown in Fig. |8](a). As expected the fidelity 
decreases when values with larger indexes are taken into 
account, since the step-like phase holograms performance 
is worse for higher-order modes. Nevertheless, values 
of F > 0.99 indicate excellent agreement with theoret- 
ical distribution. To estimate the amount of excess cor- 
relations between modes with different indexes, leading 
to n ^ > 1 in Fig. [7J one can use any appropriate 
distance measure. We choose to use vector-norm dis- 



tance D(gi%,g^) 



a {2) 



(2) 
9th, (m,n) 



Dependence plotted in Fig. [5] (b) also shows deviation 
from expected values for higher-order modes. 

In conclusion, we have investigated spatial intensity 
correlations in partially spatially coherent thermal light 
using the technique of projective measurements origi- 
nally developed to study spatial entanglement of pho- 
ton pairs. In particular, our detection scheme is similar 
to those used in experiments on orbital angular momen- 
tum entanglement (27j and essentially the same as was 
used by authors to study approximate Schmidt decom- 
position for SPDC biphotons 0- We have found much 
similarity in properties of coherent mode decomposition 
for intensity correlation function of classical thermal light 
and Schmidt decomposition for biphoton amplitude. Al- 
though, the expressions for eigenmodes and eigenvalues 
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(a) (b) 
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FIG. 8. (a) Quality of experimentally obtained distribution 
of eigenvalues. Dependence of fidelity calculated according 
to expression (|18[) on the maximal order of modes taken into 
account M + N. (b) Similar dependence of vector-norm dis- 

(2) 

tance between experimental and theoretical values of g^J n ^ , 
characterizing cross-correlations in higher-order modes. 



arc different, they are still governed by a single exper- 
imental parameter (3 = a^/ai - the ratio of coherence 
radius to the width of intensity distribution. A similar 
quantity for entangled photons is known as the Fedorov 
ratio R = Ax c /Ax s - the ratio of widths of coincidence 
(conditional) and single counts (single particle) distribu- 
tions in either coordinate or momentum space of pho- 
tons [28| . This quantity is shown to be equal to Schmidt 
number K, which is the effective number of non-zero 



eigenvalues in the Schmidt decomposition, defined as 
K = A7 2 - Thus for a pure two-photon state it can 
serve as a measure of spatial entanglement. It is clear, 
that this quantity has a purely classical analogue - opti- 
cal etendue @. Results of this work show, that analogy 
between spatial entanglement and classical spatial cor- 
relations goes further. In purely operational terms the 
qualitative similarity is almost complete, with the differ- 
ence being only quantitative - excess of g( 2 > over unity 
is limited to 1, while it can be unlimited in the case of 
entangled photons. 

One can consider the described experiments as a spe- 
cial case of thermal light ghost interference/imaging with 
mode-sensitive detection. In this context our work adds 
new arguments to the discussion of quantum and classical 
features of spatial correlations underlying this method. 
The presented results may contribute to deeper under- 
standing of what is classical and what is really quantum 
in spatial correlations of multi-mode light. 
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